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Abstract 

We present a model-independent study of spectator effects, which are 
responsible for the lifetime differences between beauty hadrons. These 
effects can be parametrized in terms of hadronic matrix elements of four 
four-quark operators. For B mesons, the coefficients of the non-factorizable 
operators turn out to be much larger than those of the factorizable ones, 
limiting considerably the usefulness of the vacuum insertion approxima- 
tion. Non-factorizable contributions to the lifetime ratio t(B~)/t(B c i) 
could naturally be of order 10-20%, and not even the sign of these contribu- 
tions can be predicted at present. In the case of the baryon, heavy-quark 
symmetry is used to reduce the number of independent matrix elements 
from four to two. In order to explain the large deviation from unity in 
the experimental result for t(Ai ) )/t(B ( i), it is necessary that these baryon 
matrix elements be much larger than those estimated in quark models. 
We have also reexamined the theoretical predictions for the semileptonic 
branching ratio of B mesons and charm counting, finding that, given the 
present theoretical and experimental uncertainties, there is no significant 
discrepancy with experiment. 
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1 Introduction 



In this paper we study "spectator effects" in inclusive decays of beauty hadrons. 
These effects involve the participation of the light constituents in the decay and 
thus contribute to the differences in the decay widths and lifetimes of different 
species of beauty hadrons. Indeed, one of our goals is to understand theoretically 
the experimental results for the lifetime ratios M: 



■{Br 



r(B d ) 

r_m 

r(B d ) 
t(A») 
r(B d ) 



1.02 ±0.04, 
1.01 ±0.07, 

0.78 ±0.05. (1) 



Here t(B s ) refers to the average 5 s -meson lifetime. Our study is performed in 
the framework of the heavy-quark expansion, in which these ratios are computed 
as series in inverse powers of the mass of the b quark (for recent reviews, 

see refs. || |7|]). The leading term of this expansion corresponds to the decay 
of a free b quark. This term is universal, contributing equally to the lifetimes 
of all beauty hadrons. Remarkably, the first correction to this result is of order 
(Aqcd/w6) 2 0, IJ- This leads to the theoretical predictions (see section § below) 



t(B- 



l + 0(l/mg), 
(1.00 ± 0.01) + 0(l/m 3 b ) 



r(B d ) 
r{B s , 



r(B d ) 



0.98 ± 0(l/m£) . (2) 



The deviation from this expectation for r(Afe)/r(i?rf) is striking, and is the prin- 
ciple motivation for this study. 

Spectator effects, i.e. contributions from decays in which a light constituent 
quark also participates in the weak process, have first been considered in refs. ||- 
p7Qf| . For decays of heavy particles, these effects are strongly suppressed due to the 
need for the b quark and a light quark in the heavy hadron to be close together (i.e. 
from a factor of the "wave-function at the origin"). As the portion of the volume 
that the b quark occupies inside the hadron is of order (Aq CD /mb) 3 , spectator 
effects appear only at third order in the heavy-quark expansion, and it might 
seem safe to neglect them altogether. However, as a result of the difference in the 
phase-space for 2 — > 2-body reactions as compared to 1 — ► 3-body decays, these 
effects are enhanced by a factor of order 16n 2 . It is conceivable that they could be 
larger than the terms of order (AqcD/ m fe) 2 included in (Q). Moreover, spectator 
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effects explicitly differentiate between different species of beauty hadrons. In 
order to understand the structure of lifetime differences, it is therefore important 
to reconsider the analysis of such effects. The striking experimental result for the 
short At, lifetime gives an additional motivation to such a study. 

Previous studies of spectator effects in the decays of beauty hadrons were 
performed using the formalism developed in refs. P|-|10|| and made two simpli- 
fying assumptions: first, the hadronic matrix elements of these operators were 



estimated employing the vacuum insertion approximation [ I I] for mesons and 
quark models || |12| for baryons, and secondly the mass of the charm quark was 
neglected in the calculation of the coefficients of the four-quark operators in the 
heavy-quark expansion. In order to explore the acceptable range of theoretical 
predictions, we do not impose factorization or quark-model approximations on 
the hadronic matrix elements. Instead, we parametrize them by a set of hadronic 
parameters and see how the lifetime ratios depend upon them. One of our main 
conclusions is that only a detailed field-theoretic calculation of the relevant ma- 
trix elements can lead to reliable predictions. In addition, we derive the exact 
expressions for the coefficients as functions of the charm-quark mass. 

The semileptonic branching ratio of B mesons has also received considerable 
attention. For many years it appeared that the theoretical predictions for this 
quantity jl3|-||l5| lay above the measured value. More recently, the theoretical 



predictions have been refined by including exact expressions for the 0(a s ) cor- 
rections fll6| . Using the results of these calculations, we present a new analysis 
of the semileptonic branching ratio Ssl and the average charm multiplicity n c in 
B decays. We find that the freedom in the choice of the renormalization scale 
(which reflects the ignorance of higher-order perturbative corrections) allows us 
to obtain consistent predictions for both quantities simultaneously. We also cal- 
culate the spectator contributions to B^h and n c and show that they could change 
the semileptonic branching ratio by an amount of order 1%, whereas their effect 
on n c is negligible. 

In section we discuss the heavy-quark expansion for inclusive decay rates, 
and we present our results for the contributions arising from spectator effects. 
In section |3|, we introduce a set of hadronic parameters defined in terms of the 
relevant four-quark operator matrix elements between 5-meson and A^-baryon 
states. Heavy-quark symmetry is used to derive some new relations between 
the baryonic matrix elements of these operators. In section ^, we then discuss 
the phenomenological implications of our results for the understanding of beauty 
lifetimes. We also present a critical discussion of previous estimates of spectator 
effects based on the factorization approximation. A detailed discussion of the 
semileptonic branching ratio of B mesons is presented in section |5]. Section ^| 
contains the conclusions. The renormalization of the operators and parame- 
ters describing the spectator effects, and the behaviour of these parameters with 
respect to the large-A^ c limit, are discussed in appendix A, while appendix B 
contains details about the calculation of the semileptonic branching ratio. 
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The reader who is primarily interested in the phenomenological implications 
of our analysis can omit sections |] and |3] and proceed directly to sections f|, [| 
and |6|, which are written in a self-contained way. 



2 Heavy- quark expansion 

Inclusive decay rates, which determine the probability of the decay of a particle 
into the sum of all possible final states with a given set of quantum numbers 
{/}, have two advantages from the theoretical point of view: first, bound-state 
effects related to the initial state can be accounted for in a systematic way using 
the heavy-quark expansion; secondly, the fact that the final state consists of a 
sum over many hadronic channels eliminates bound-state effects related to the 
properties of individual hadrons. This second feature is based on the hypothesis 
of quark-hadron duality, i.e. the assumption that cross sections and decay rates 
are calculable in QCD after a "smearing" procedure has been applied [[HJ. We 



shall not discuss this hypothesis here; however, if after the non-perturbative eval- 
uation of the spectator effects discussed in our analysis there remain significant 
discrepancies between theory and experiment (for the lifetime ratio T(Ab)/r(Bd), 
in particular), one may have to seriously question the assumption of duality. A 
recent study of inclusive B decays, in which duality violations are invoked to add 
non-perturbative contributions of order Aqcd / not present in the heavy-quark 
expansion, can be found in ref. [IE\ . 



Using the optical theorem, the inclusive decay width of a hadron H b containing 
a b quark can be written as the forward matrix element of the imaginary part of 
the transition operator T, 

T(H b ^X) = —Im (H b \ T \H b ) = — !— (H b \ T \H b ) , (3) 
m H , 2m Hh 



where T is given by 



d 4 xT{£ eff (x),/: efr (0)}. (4) 



For the case of semileptonic and non-leptonic decays, the effective weak La- 
grangian, renormalized at the scale \i = m b , is 

£ c s = -i=- V ch \c x {m h ) d' L ^^u L ci/fb L + s'^^cl cj/fb L 



+ c 2 (m b ) [cl^^ul d' L Yb L + c L ^^c L s' L ^b L 

+ E hj^i CL^bA + h.c, (5) 

0—n ,, n- > 
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where qi = §(1 — 7s)<? denotes a left-handed quark field, d! = d cos9 c + s sin# c 
and s f = s cos6 c — d sin6 c are the Cabibbo-rotated down- and strange-quark 
fields (sin# c ~ 0.2205), and we have neglected b —> u transitions. The Wilson 
coefficients c\ and C2 take into account the QCD corrections arising from the fact 
that the effective Lagrangian is written at a renormalization scale /i = m b rather 
than mw They can be calculated in perturbation theory. The combinations 
c± = C\ ± C2 have a multiplicative evolution under change of the renormalization 



scale. To leading order, they are given by [ 19]— [21 



fa s {m w )\ 12 
c±(m 6 ) = — - — - , a_ = -2a + = -— — — . 6 
\ a s (m b ) J 33 - 2n f 

In the numerical analysis we shall take the values c + (rrib) — 0.86 and C-(m b ) = 
l/c+(m b ) ~ 1.35, corresponding to a s (mz) = 0.117. 

Since the energy release in the decay of a b quark is large, it is possible 
to construct an Operator Product Expansion (OPE) for the bilocal transition 
operator (|j), in which it is expanded as a series of local operators with increasing 
dimension, whose coefficients contain inverse powers of the 6-quark mass. The 
operator with the lowest dimension is bb. There is no independent operator with 
dimension four, since the only candidate, bipb, can be reduced to bb by using 
the equations of motion ||, |J. The first new operator is bg^^G^b and has 
dimension five. Thus, any inclusive decay rate of a hadron Hf, can be written in 
the form 

vm Y . G\m\ 1 f f,„ ]lh] „. , f (H b \bg s a, u G^b\H b ) \ 

(7) 

where c{ are calculable coefficient functions (which also contain the relevant CKM 
matrix elements) depending on the quantum numbers {/} of the final state. For 
semileptonic and non-leptonic decays, the coefficients c{ have been calculated at 
one-loop order |22|, |23], |16| , and the coefficients c( at tree level || [24|. 



In the next step, the forward matrix elements of the local operators in the 
OPE are systematically expanded in inverse powers of the 6-quark mass, using 
the heavy-quark effective theory (HQET) pSJ . One finds & H| 



2m 



— (H b \ bb \H b ) - 1 ^ + 0(l/m b ) , 

1 (H b \ bg.<T IB ,G> a, b \H b ) = 2fi 2 G (H b ) + 0(l/m b ) , (8) 



where ^{H b ) and fiQ(H b ) parametrize the matrix elements of the kinetic-energy 
and the chromo-magnetic operators, respectively. The purpose of doing this ex- 
pansion is that whereas the matrix elements in (0) contain an implicit dependence 
on the 6-quark mass, the parameters appearing on the right-hand side of (M) are 
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independent of m& (modulo logarithms). These parameters can be determined, 
to some extent, from the spectrum of heavy hadron states. Below we shall need 
the values 

/4(A 6 ) - (4(B) = -(0.01 ± 0.03) GeV 2 , 

/4(5) = ~(m**-m|)~0.36 GeV 2 , 
f x 2 G (A b ) = 0. (9) 
The difference /x 2 (A b ) — (4(B) can be extracted from the mass formula [|. [?j 



{m Ab - m Ac ) - [m B - m D ) 



/4(B)-/4(A.)](^-^) + O(l /n 4),(10) 



where m B = /(wb + 3ms*) and = 7 (mc + 3hid.) denote the spin- averaged 



meson masses. With m Ab = (5625 ± 6) MeV f26]j , this relation leads to the value 
quoted above. 

To order 1/m 2 in the heavy-quark expansion, the lifetime ratio for two beauty 
hadrons is given by 

7(^f) + 2^ + CG ml + (11) 

where cq ~ 1.2 can be obtained using the results of refs. || 24j . Using then 



the values given in (|9]), and assuming that in the case of the B s meson SU(3)- 
breaking effects in the values of the matrix elements are of order 20%, we arrive at 
the predictions given in (fj). Note that in taking a ratio of lifetimes, theoretical 
uncertainties related to the values of the 6-quark mass (including renormalon 
ambiguities) and CKM elements cancel to a large extent. It is for this reason 
that we restrict our discussion to the calculation of ratios of lifetimes and decay 
rates. 

The first two terms in the heavy-quark expansion (|?|) arise from decays in 
which the (light) spectator quarks interact only softly. Additional contributions 
of this type appear at the next order through gluonic operators of dimension 
six, such as b , y f j - (iD u G tJ-u )b. Since matrix elements of these operators are blind 
to the flavour of the spectator quarks, they can be safely neglected in our anal- 
ysis. "Hard" spectator effects manifest themselves first in the matrix elements 
of four-quark operators of dimension six. Some examples of the corresponding 
contributions to the transition operator T are shown in figure |l[ They only 
appear in the heavy-quark expansion of non-leptonic decay rates.Q Since these 
contributions arise from one-loop rather than two-loop diagrams, they receive a 
phase-space enhancement factor of order 167r 2 relative to the other terms in the 
heavy-quark expansion. 

lr This is no longer true if b — > u transitions or the decays of the B c meson are considered. 
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bTjqq Ti b 




Figure 1: Spectator contributions to the transition operator T (left), and the 
corresponding operator in the OPE (right). Here Tj denotes some combination 
of Dirac and colour matrices. 



We have calculated the coefficients of the corresponding four-quark operators 
at tree level, including for the first time the dependence on the mass of the 
charm quark. This extends the results obtained in ref. 0. We find that the 
corresponding contributions to the non-leptonic widths of mesons and baryons 
containing a h quark are given by the matrix elements of the local operator 



spec 



2G 2 F m 2 b 

71 

IGlml 



cb\ 



:i-z) 



2 ClC2 + — (c? + 4) ) oi_ A + 2(4 + 4) t;-_ a 



3tt 



i + 



z 



O d ;_ A -(l + 2z)O d 



S-P 



+ 24 



l + -)n_ A -(l + 2z)Tl_ P 



2G 2 F m 2 h , l9 , if 1 9 9 

\V cb \ 2 I [2c lC2 + — 4 + N.4 



+ 24 



'l-z)T, 



V-A 



l-z)0 
1 + 2z) Tg_ P 



v _ A -(l + 2z)O s s _ P 
(12) 



where z = m^/m 2 , and iV c =3 is the number of colours. The local four-quark 
operators appearing in this expression are defined by 

O q v _ A = hj^qL q L l^b L , 
O q s _ P = b R q L q L b R , 



Ty-A = b L l^t a q L q L l^t a b L 
T|_ p = b R t a q L q L t a b R , 



(13) 



where t a = \ a /2 are the generators of colour SU(3). For dimensional reasons, 
r spcc is proportional to m 2 rather than mf, in accordance with the fact that 
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spectator effects contribute at third order in the heavy-quark expansion. The 
first term in (|I2|) arises from the upper diagram in figure p], whereas the second 
and third terms come from the contributions of the lower diagram with a cu and 
cc quark pair in the loop. We note that in the limit z = our results agree with 
ref. 0, and with the corresponding expression derived for the lifetimes of charm 
hadrons in refs. [0-111011. 



The operators in ([13]) are renormalized at the scale m&, which will be im- 
plicit in our discussion below. This choice has the advantage that logarithms 
of the type [a s ln(m&///had)] n , where /ihad is a typical hadronic scale, reside en- 
tirely in the hadronic matrix elements of the renormalized operators. Using the 
renormalization-group equations, the expressions presented in this paper can be 
rewritten in terms of operators renormalized at any other scale. However, at 
present the scale dependence of the renormalized operators below the scale m;, is 
known only to leading logarithmic order [E7], ESJ. It is discussed in detail in ap- 



pendix A. Since here we shall treat the matrix elements as unknown parameters, 
we can avoid all uncertainties related to the operator evolution by working at the 
scale m b . 

The hadronic matrix elements of the four-quark operators in contain 
the non-perturbative physics of the spectator contributions to inclusive decays of 
beauty hadrons. However, the same operators also contribute to the decay of the 
b quark, through tadpole diagrams in which the light-quark fields are contracted 
in a loop. These "non-spectator" contributions are independent of the flavour of 
the light quark q and thus contribute equally to the decay widths of all beauty 
hadrons. They are not of interest to our discussion here. In order to isolate the 
true spectator effects, we shall implicitly assume a normal ordering of the four- 
quark operators, which has the effect of subtracting tadpole-like diagrams. In 
practice, this is equivalent to choosing a particular renormalization prescription 
for the operators. Alternatively, one may isolate the spectator effects by con- 
sidering light-quark flavour non-singlet combinations of the operators; thus, for 
example, for matrix elements between states one could take (O d - O u ) and 
{T d - T u ) instead of O d and T d . 



3 Parametrization of the matrix elements 

In previous analyses ||, [Pl-jlOfl, the hadronic matrix elements of the four-quark 



operators in (|T3|) have been estimated making simplifying assumptions. Here 
we shall avoid such assumptions and, without any loss of generality, express the 
relevant matrix elements in terms of a set of hadronic parameters. Clearly, such 
an approach has less predictive power; however, it does allow us to find the 
range of predictions that can be made in a model-independent way. In view of 
the apparent discrepancy between theory and experiment for the Af, lifetime, we 
find it worth while to question the model-dependent assumptions made in earlier 



S 



analyses. Ultimately, the relevant hadronic parameters may be calculated using 
some field-theoretic approach such as lattice gauge theory or QCD sum rules. It 
has also been suggested that combinations of these parameters may be extracted 
from a precise measurement of the lepton spectrum in the endpoint region of 
semileptonic B decays, or from a study of spectator effects in charm decays p| . 



2m Bo 



2m Bc 



2m Ba 



3.1 Mesonic matrix elements 

For matrix elements of the four-quark operators between 5-meson states, we 
define parameters Bi and Si such that: 

' / R I r)Q I R \ — ^ B q mB i r> 
\ri q \ Uy_ A \Bg) = B x , 

(B q \Ol_ P \B q )^ f -i^B 2 , 

(B q \n_ A \B q )^ f -i^e 1 , 

^(B q \T Lp \B q ) = ^e 2 . (14) 

This definition is inspired by the vacuum insertion (or factorization) approxi- 
mation ||11|| , according to which the matrix elements of four-quark operators are 
evaluated by inserting the vacuum inside the current products. This leads to 

(B,\ Ol_ A |B t ) = r^l) 2 <B„| 0|_ P \B q ) = , 
(B q \T*_ A \B q ) = (B q \T$_ P \B q )=0, (15) 

where fa is the decay constant of the B q meson, defined as 

(0\qn 5 b\B q (p))=tf BqP r (16) 

Hence, the factorization approximation corresponds to setting Bi = 1 and = 
at some scale \x (which in general will be different from our adopted choice fi = 
rrib), where the approximation is believed to be valid. The exact values of the 
hadronic parameters are not yet known. However, as discussed in appendix A, 
in the large- N c limit 

B % = 0(1) , = 0(1/N C ) . (17) 

An estimate of the parameters £j using QCD sum rules has been obtained by 



Chernyak, who finds that e\ ~ —0.15 and e 2 ~ [30 
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In terms of the parameters Bi and £j, the matrix elements of the operator 
r spec in dH) are: 

— (B-\ r spcc \B~) = r r/ S pcc (1 - zf {(2c 2 , - c 2 _) Ex + 3(4 + c 2 ) ex} , 



2m b 
1 



2m 



(■Srfl r spec |-Bd) 



1 



-r ^spcc (1 - ^) 2 cos 2 6> c ^ - (2c + - c_) 2 



1 



1 + - ) B x - (1 + 2z) B 2 



+ - (c+ + c. 



l + -)£i-(l + 2*)e2 



- T r/ spcc Vl^4i sin 2 # c j - (2c + - c_) 2 [(1 - z) B x - (1 + 2z) B 2 



+ - (c+ + c. 



where c± = ci ± c 2 , and 



'l-z)ei-(l + 2z)e 2 



G 2 F ml 

1927T 3 



77 spec = 167T' 



ni) 



(19) 



The spectator contribution to the width of the B s meson is obtained from that 
of the B d meson by the replacements sin^c <-> cos^ c and fB, m B fB s i m B s - Of 
course, the values of the parameters Bi and Si for the B s meson will also differ 
from those for the Bd meson due to SU(3)-breaking effects. 

Two remarks are in order regarding the result dl8|) . The first concerns the 
expected order of magnitude of the spectator contributions to the total decay 
width of a B meson. At leading order in the heavy-quark expansion, the total 
width of a beauty hadron is r tot ~ 3.7 x r , where the numerical factor arises 
from the phase-space contributions of the semileptonic and non-leptonic channels 
0.085) 



v we use z 



2J. It follows that 



I\ 



Vs 



tot 



5%. 



(20) 



The second remark concerns the structure of the coefficients in fll8|) . Given 
that c + ~ 0.86 and c_ ~ 1.35, one observes that the coefficients of the colour 
singlet-singlet operators are one to two orders of magnitude smaller than those 
of the colour octet-octet operators. This implies that at the scale m& even small 
deviations from the factorization approximation can have a sizeable impact on 
the results. 
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3.2 Baryonic matrix elements 

Next we study the matrix elements of the four-quark operators between A;,-baryon 
states. Since the is an iso-singlet, the matrix elements of the operators with 
q = u or d are the same, and below we drop this label. In the case of baryons, 
we find it convenient to use the colour identity 

(t a )a(3 (t a ) y S = - S a s 5 7/ 3 — 5 a f3 <5 7< 5 (21) 

to rewrite T = — ~ O + ^O and introduce the operators are colour indices) 
Ov-A = \l^ L &fV L , Os-p = ¥ R q{ qi b R . (22) 

instead of Ty-A and Ts-p. 

The heavy-quark spin symmetry, i.e. the fact that interactions with the spin 
of the heavy quark decouple as the heavy-quark mass tends to infinity, allows us 
to derive two relations between the matrix elements of the four-quark operators 
between A^-baryon states. To find these relations, we note that the following 
matrix element vanishes in the limit m;, — > oo: 

(Afc| BVtt V q k L l"q l L |A*> = 0(l/m b ) . (23) 



2m A . 



The physical argument for this is that, because of the spin symmetry for heavy 
quarks, the matrix elements for left-handed and right-handed b quarks must be 
the same. The above result then follows since &7 M 7s& = bp/y^R — 6l7/x&l- A 
more formal argument can be given using the covariant tensor formalism of the 
HQET p5l [IT] to show that the matrix element in ( p3|) is proportional to 



(2S& • Slight) = S\ b (SA b + 1) — Sb(Sb + 1) — Slight (-Slight + 1) = , (24) 

since, in the heavy-quark limit, the light degrees of freedom are in a state with 
total spin zero. 

Using the Fierz identity 

^7m75 V q k L ^q[ = -2 b R q[ q\ V R - U Ll ,q l L q k L YK , (25) 
we then obtain the relations 

' (A 6 | Os-p |A 6 ) = -]- (A 6 | Ov^a |A 6 ) + 0{l/m b ) , 



2m Ab 2 2m Ab 

(A 6 | Os-p |A 6 ) = -\ (A 6 ] Ov-A \A b ) + 0(l/m b ) . (26) 
2m\ b Z zm Ab 

The corrections of order l/m b to these relations contribute at order 1 / in the 
heavy-quark expansion and so are negligible to the order we work in. This leaves 
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us with two independent matrix elements of the operators Oy-A and Ov-A- The 
analogue of the factorization approximation in the case of baryons is the valence- 
quark assumption, in which the colour of the quark fields in the operators is 
identified with the colour of the quarks inside the baryon. Since the colour wave 
function for a baryon is totally antisymmetric, the matrix elements of Oy-A and 
Oy-A differ in this approximation only by a sign. Hence, we define a parameter 
B by 

(A 6 | 6 V _ A \A b ) = -B (A*] O v . A |A 6 ) , (27) 

with B = 1 in the valence-quark approximation. 

For the baryon matrix element of Ov-a itself, our parametrization is guided 
by the quark model. We write 

(A b \0 V - A \A b ) = -3^r, (28) 



2m Ab n 1 48 

where in the quark model r is the ratio of the squares of the wave functions 
determining the probability to find a light quark at the location of the b quark 
inside the A b baryon and the B meson, i.e. |], I(| 



W(Q)I 
l< 9 (o)| 



- = ^e-^- (29) 



Guberina et al. have estimated the ratio r for charm decays and find r ~ 0.2 in 
the bag model, and r ~ 0.5 in the non-relativistic quark model ||. This latter 
estimate has been used in more recent work on beauty decays ||. A similar 
result, r ~ 0.1-0.3, has been obtained by Colangelo and Fazio using QCD sum 
rules [[32]. Recently, Rosner has criticized the existing quark- model estimates of 
r. Assuming that the wave functions of the A& and Sj baryons are the same, 
he argues that the wave-function ratio should be estimated from the ratio of the 



spin splittings between E fe and baryons and B and B* mesons p3| . This leads 
to 

r = \^ ^. (30) 

3 m z B , — m z B 



If the baryon splitting is taken to be m|, — m| t ~ m|* — m| c = (0.384 ± 

0.035) GeV 2 , this leads to r ~ 0.9 ± 0.1. If, on the other hand, one uses the 
preliminary result m^* — ms t = (56 ± 16) MeV reported by the DELPHI Collab- 



oration f34| , one obtains r ~ 1.8 ± 0.5. We conclude that it is conceivable that 
r ~ 1, i.e. larger than previous estimates. 

In terms of these parameters, the matrix element of r spec is given by 



1 r f 

(A 6 | r spcc |A 6 ) = r 77 spec — l 4(1 - z) 



r\ \--u\ spec u " /spec 

2m A , 16 



[c 2 _ - 4) + (c 2 _ + 4) B 



12 



1 - zf(l + z) cos 2 9 c + y/l-Az sin 9, 
(c_ - c+) (5c+ - c_) + (c_ + c+) 2 B 



(31) 



3.3 Numerical results 

To illustrate the main features of our results, we calculate the coefficients of 
the hadronic parameters Bi, Si, and r and Br in the matrix elements fll8|) and 
fl3"ID in units of To^spec- In order to study the dependence on the mass ratio 
z = (m c /mb) 2 , we first keep the values of the Wilson coefficients in the effective 
Lagrangian fixed and vary the mass ratio in the range z = 0.085 ± 0.015. This 
leads to the numbers shown in table [TJ, where the variation with z is indicated 
as a change in the last digit (s). Note that for mesons the coefficients of the 
parameters Bi are much smaller than those of the parameters £j. It is apparent 
that the results are rather stable with respect to the precise value of z. From 
now on we shall always use the central value z = 0.085, which is obtained, for 
instance, for m c = 1.4 GeV and rn& = 4.8 GeV. 



Table 1: Coefficients of the hadronic parameters obtained for z = 0.085 ± 
0.015. The values c + = 0.861 and c_ = 1.349 are kept fixed. 





B 1 


B 2 


£i 


£2 


r 


Br 


B~ 


-0.28(1) 




6.43(21) 








B d 


-0.04(0) 


0.05(0) 


-2.12(6) 


2.39(2) 






B s 


-0.03(0) 


0.04(0) 


-1.83(11) 


2.32(5) 






A b 










0.14(1) 


0.26(1) 



As another check on the stability of the results we present in table |2| the 
coefficients obtained using different "matching" procedures. To this end, we 
renormalize the Wilson coefficients c+ and c_ of the effective Lagrangian at a 
scale fi different from m&. Thus c + and c_ are modified by replacing mj by // 
in @. This gives us predictions in terms of operators renormalized at /i, which 
we then rewrite in terms of those defined at m& by using the evolution equations 
given in appendix A. Thus, the meaning of the hadronic parameters is the same 
as before, and the differences between the numerical results can be viewed as 
an estimate of unknown higher-order perturbative corrections, which we neglect 
throughout this paper. The coefficients of the parameters Bi for mesons, as 
well as of the parameters r and Br for the A& baryon, show a significant scale 
dependence. To reduce this dependence would require a full next-to-leading order 
calculation of radiative corrections, which is beyond the scope of this paper. 
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Table 2: Coefficients of the hadronic parameters obtained for the matching 
scales \i = m b /2, m b and 2m b , with m b = 4.8 GeV. The value z = 0.085 is kept 
fixed. 







B\ B 2 Ei e 2 


r Br 


B~ 


m b j l 
m b 
2m b 


— u.oz — (.ZD — 
-0.28 6.43 
+0.02 5.90 




B d 


m b /2 
m b 
2m b 


-0.04 0.04 -2.32 2.61 
-0.04 0.05 -2.12 2.39 
-0.08 0.09 -1.98 2.24 




B s 


m b /2 
m b 
2m b 


-0.03 0.04 -2.00 2.54 
-0.03 0.04 -1.83 2.32 
-0.07 0.08 -1.72 2.18 




A b 


m b /2 
m b 
2m b 




0.21 0.30 
0.14 0.26 
0.09 0.23 



4 Phenomenology of beauty lifetimes 

We shall now discuss the phenomenological implications of our results for the 
calculation of beauty lifetime ratios. The spectator contributions to the decay 
widths of B q mesons and of the A b baryon are described by a set of hadronic 
parameters: Bi 2 and for mesons, and r and B for baryons. The explicit 
dependence of the decay rates on these quantities is shown in and (pip. 
The numerical values of the coefficients multiplying the hadronic parameters are 
given in table for three different choices of the matching scale /x. For the 



numerical analysis we need the value of the parameter r] spec defined in (|19|) . We 
take fs = 200 MeV and m b = 4.8 GeV, so that 77 spe c — 0.30, and absorb the 
uncertainty in rj spec into the values of the hadronic parameters. 



4.1 Lifetime ratio r(B )/r(B d ) 

We start by discussing the lifetime ratio of the charged and neutral B mesons. 
Because of isospin symmetry, the lifetimes of these states are the same at order 
l/m b in the heavy-quark expansion, and differences arise only from spectator 
effects. If we write 

T -±%\ = 1 + fciSi + k 2 B 2 + k 3£l + he 2 , (32) 

the coefficients ki take the values shown in table |3|. The most striking feature of 
this result is the large imbalance between the coefficients of the parameters Bi 
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and £j, which parametrize the matrix elements of colour singlet-singlet and colour 
octet-octet operators, respectively. With of order l/N c , it is conceivable that 
the non-factorizable contributions actually dominate the result. Thus, without a 
detailed calculation of the parameters Si no reliable prediction can be obtained. 
In this conclusion we disagree with the authors of ref. ||, who use factorization 
(at a low hadronic scale) to argue that r(B~) /t{B^) must exceed unity by an 
amount of order 5%. We will return to this in section f4~4| below. 



Table 3: Coefficients ki appearing in 



n 


h k 2 h h 


m b /2 
m b 
2rrib 


+0.044 0.003 -0.735 0.201 
+0.020 0.004 -0.697 0.195 
-0.008 0.007 -0.665 0.189 



The experimental value of the lifetime ratio given in ([I]) can be employed to 
constrain a certain combination of the parameters: 



£l = T 



(0.02 ± 0.04) - k 1 B 1 - k 2 B 2 - k 4 e 2 ] . (33) 

With the values of the coefficients given in table |3], this relation implies 

£i~0.3e 2 + (5, (34) 

where we expect \S\ < 0.1. This becomes a particularly useful constraint if the 
parameters £j turn out to be large. Below, we shall use ( |33"D to eliminate e\ from 
our predictions for spectator effects. 



4.2 Lifetime ratio r(B s )/r(Bd) 

In the limit where SU(3)-breaking effects are neglected, the spectator contribu- 
tions to the decay widths of B s and Bj mesons are too similar to produce an 
observable lifetime difference (see table |2|). The corresponding contributions to 
the lifetime ratio t{B s )/t{B ( j) are of order 10~ 3 . A precise prediction for SU(3)- 
breaking effects is difficult to obtain. Allowing for 30% SU(3)-breaking effects 
in the matrix elements of the four-quark operators describing the spectator con- 
tributions, we estimate that the resulting contributions to the lifetime ratio are 
of order 1-2%. Contributions of order 1% (or less) could also arise from SU(3)- 
breaking effects in the matrix elements appearing at order 1 / in the expansion 
(ITT]). Hence, we agree with ref. that 



t(B.) 
r(B d ) 



1 ± 0(1 ( , 



(35) 
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4.3 Lifetime ratio r(A b )/r(B d ) 

As mentioned in the introduction, the low experimental value of the lifetime ratio 
r{h.b)/r{Bd) is the primary motivation for our study. We shall now discuss the 
structure of spectator contributions to this ratio. It is important that heavy-quark 
symmetry allows us to reduce the number of hadronic parameters contributing 
to the decay rate of the A& baryon from four to two, and that these parameters 
are almost certainly positive (unless the quark model is completely misleading) 
and enter the decay rate with the same sign. Thus, unlike in the meson case, the 
structure of the spectator contributions to the width of the A& baryon is rather 
simple, and at least the sign of the effects can be predicted reliably. 

For a more detailed discussion, we distinguish between the two cases where 
one does or does not allow spectator contributions to enhance the theoretical 
prediction for the semileptonic branching ratio, -Bsl, of B mesons. As we will 
discuss below, the theoretical prediction for .Bsl, which neglects spectator con- 
tributions, is slightly larger than the central experimental value. If spectator 
effects increased the prediction for .Bsl further, this discrepancy could become 
uncomfortably large. 

If we do not allow for an increase in the value of the semileptonic branching 
ratio, the explanation of the low value of T{h.\ } )/r{Bd) must reside entirely in a 
low value of the A& lifetime (rather than a large value of the 5-meson lifetime). 
This can be seen by writing 

r(A b ) _ , A , (r(B-)\ 1/2 1 r(A 6 ) ( r( B -)\ 1/2 

(36) 

where -Bsl is the average semileptonic branching ratio of B mesons, and Fsl(-B) 
is the semileptonic width. In the second step we have replaced the geometric 
mean [t(B~) r(B d )] 1/2 by the average .B-meson lifetime, which because of isospin 
symmetry is correct to order l/mf in the heavy-quark expansion. Since there 
are no spectator contributions to the semileptonic rate Fsi(B), and since we 
do not allow an enhancement of the semileptonic branching ratio, in order to 
obtain a small value for t(A&)/t(.BA we can increase the width of the A& baryon 
and/or decrease (within the experimental errors) the lifetime ratio t(B~)/t(B ( i). 
Allowing for a downward fluctuation of this ratio by two standard deviations, 
i.e. T(B~)/r(Bd) > 0.94, and using the estimate of l/m| corrections in (g), we 
conclude that 



IM > 0.97 x 



0.98 - rs P ec ( Afe ) 



T(A 6 



0.95- (d 1 + d 2 B)r, (37) 



where r spec (A(,) is the spectator contribution to the width of the A& baryon. The 
values of the coefficients di are given in table |]. If we assume that r and B 
are of order unity, we find that the spectator contributions yield a reduction of 
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the lifetime of the A& baryon by a few per cent, and that t(Ai,)/t(B ( i) > 0.9, 
in contrast with the experimental result given in ([I]). If, for example, we try to 
push the theoretical prediction by taking the large value B — 1.5 (corresponding 
to a violation of the valence-quark approximation by 50%) and choosing a low 
matching scale fi = m b j2, we have to require that r > r min with r min = 3.1, 
2.2 and 1.3 for t(Ai ) )/t(B ( j) = 0.78, 0.83 and 0.88 (corresponding to the central 
experimental value and the la and 2a fluctuations). Hence, even if we allow for 
an upward fluctuation of the experimental result by two standard deviations, we 
need a value of r that is significantly larger than most quark-model predictions 
(see the discussion in section |3.2|) . Clearly, a reliable field-theoretic calculation 
of the parameters r and B is of great importance to support or rule out such a 
possibility. 



Table 4: Coefficients di appearing in (|37D and (j38|). 





d\ d 2 d$ G?4 


m b /2 
m b 
2m b 


0.016 0.023 0.178 -0.201 
0.012 0.021 0.173 -0.195 
0.008 0.020 0.167 -0.189 



On the other hand, the low experimental value of the semileptonic branching 
ratio may find its explanation in a low renormalization scale (see section [| below), 
or it may be caused by the effects of New Physics, such as an enhanced rate for 



flavour- changing neutral currents of the type b — > sg |p5|— 1|38|| . Hence, one may 
be misled in using the semileptonic branching ratio as a constraint on the size 
of spectator contributions. Then there is the possibility to decrease the value of 
r(Afe)/r(i? ( i) by increasing the lifetime of the B^ meson, i.e. in fl37|) we can allow 
for spectator contributions to the width of the B^ meson. From table ^ it follows 
that the contributions of the parameters B\ and B 2 are very small (of order 10~ 3 ) 
and can safely be neglected. Thus, we obtain 

4^4 ^ 0.98 - (d 1 + d 2 B) r - (tfeei + d 4 e 2 ) , (38) 

where the coefficients d% and d 4 are also shown in table (|. At first sight, it 
seems that with a positive E\ and a negative e 2 of order l/N c one could gain 
a contribution of about —0.1, which would take away much of the discrepancy 
between theory and experiment. However, the experimental result for the lifetime 
ratio t(B~)/t(B,i) imposes a useful constraint. Using (|33|) to eliminate E\ from 
the relation (|38l), and allowing the parameters Bi to take values between and 
2, we find 

444 - °- 98 ± °- 02 + °- 15£ 2 - {di + d 2 B) r > 0.88 - (di + d 2 B) r . (39) 

T \ B d) 
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where in the last step we have assumed that \e2\ < 0.5, which we consider to be 
a very conservative bound. Even in this extreme case, a significant contribution 
must still come from the parameters r and B. 

In view of the above discussion, the short A b lifetime remains a potential prob- 
lem for the heavy-quark theory. If the current experimental value persists, there 
are two possibilities: either some hadronic matrix elements of four-quark opera- 
tors are significantly larger than naive expectations based on large- N c counting 
rules and the quark model, or (local) quark-hadron duality, which is assumed 
in the calculation of lifetimes, fails in non-leptonic inclusive decays. In the sec- 
ond case, the explanation of the puzzle lies beyond the heavy-quark expansion. 
Let us, therefore, consider the first possibility and give a numerical example for 
some possible scenarios. Assume that \i = m b /2 is an appropriate scale to use 
in the evaluation of the Wilson coefficients, and that B = 1.5. Then, to obtain 
r(Ab)/r(5(i) = 0.8 without enhancing the prediction for the semileptonic branch- 
ing ratio requires r ~ 3, i.e. several times larger than quark-model estimates. If, 
on the other hand, we consider r = 1.5 as the largest conceivable value, we need 
e 2 ~ —0.5, corresponding to a rather large matrix element of the colour-octet op- 
erator Ts_p. Such a value of e 2 leads to an enhancement of the 5-meson lifetime, 
and hence to an enhancement of the semileptonic branching ratio of B mesons, 
by A£>sl — 1%- As we will discuss in section || this is still tolerable provided yet 
unknown higher-order corrections confirm the use of a low renormalization scale. 
Although in both cases some large parameters are needed, we find it important to 
note that until reliable field-theoretic calculations of the matrix elements of four- 
quark operators become available, a conventional explanation of the A^-lifetime 
puzzle cannot be excluded. 



4.4 Relation with the conventional factorization approach 

We now discuss in detail the relation of our approach with previous analyses 
based on the factorization approximation ||. This approximation amounts to 
setting (see appendix A) 

/ / \ \ 4 /A) 

£^had) = , "\ , Si{fi had ) = , (40) 



at some hadronic scale /^ad <C mb- Here f3o is the first coefficient of the j3 
function. The evolution of the operators from m& to /^ad is done to leading 



logarithmic order [[10], |27], |28|]. For the purpose of our discussion, we fix /ihad 
such that a s (/Xhad) = 0.5. Then the relation between the hadronic parameters 
renormalized at the scale m b with those renormalized at /ihad, as given in ( |A.5| ), 
is: 

Bi(m b ) ~ 1.48-B;(// ha d) - 0.36ei(/i had ) , 
£i (m b ) ~ -0.08^Gu had ) + 1.06 £i (/i had ) . (41) 
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The theoretical prediction for the lifetime ratio t(B~)/t(B,i) in terms of the 
hadronic parameters renormalized at the scale /ihad is obtained by combining ( [41)) 
with the numbers given in table |3|. We find 

" 1 + 0.085! (/i had ) - 0.015 2 (/i had ) - 0.75 ei (/i had ) + 0.20e 2 (^ had ) . (42) 

The factorization approximation -Bj(/ihad) — 0.68 and £j(/ihad) = leads to 
t(B~)/t(Bci) ~ 1.05, which is close to the value obtained by Bigi et al. 0. 
However, it is evident from (^) that even at a low scale the coefficients of the 
non-factorizable terms are still much larger than those of the factorizable ones. 
Hence, it remains true that non-factorizable corrections are potentially important. 
It would be justified to neglect these contributions only if the parameters £j(/ihad) 
were significantly smaller than 10%. However, to the best of our knowledge there 
is currently no compelling argument to support such a strong restriction. 



5 Semileptonic branching ratio of B mesons 

The semileptonic branching ratio of B mesons has received considerable attention 
in the past. It is defined as 

B T(B^Xeu) 

SL E*r(s^x£p) + r NL + r rare ' 1 ; 

where Tnl and r rare are the inclusive rates for non-leptonic and rare decays, 
respectively. The status of the experimental results on the semileptonic branching 
ratio is controversial, as there is a discrepancy between low-energy measurements 
performed at the T(4s) resonance and high-energy measurements performed at 
the Z° resonance. The average value at low energies is Bsl = (10.37 ± 0.30)% 
whereas high-energy measurements give Bj£ = (11.11 ±0.23)% dp]. The 
superscript (b) indicates that this value refers not to the B meson, but to a 
mixture of b hadrons (approximately 40% B~ , 40% B d , 12% B s , and 8% A b ). 
Assuming that the corresponding semileptonic width is close to that of the 
B meson,| we can correct for this and find -Bsl — i T {B) /r{b)) B^ = (11.30 ± 
0.26)%, where rib) = (1.57 ± 0.03) ps is the average lifetime corresponding to 
the above mixture of b hadrons [|TJ. The discrepancy between the low-energy 
and high-energy measurements of the semileptonic branching ratio is therefore 
larger than 3 standard deviations. If we take the average and inflate the error to 
account for this fact, we obtain 

B SL = (10.90 ±0.46)%. (44) 



2 Theoretically, this is expected to be a very good approximation. 
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In understanding this result, an important aspect is charm counting, i.e. the 
measurement of the average number n c of charm hadrons produced per B decay. 
Theoretically, this quantity is given by 

n c = 1 + B(B -> X ccV ) - B(B -> no charm) , (45) 

where B(B — > X cgs /) is the branching ratio for decays into final states containing 



two charm quarks, and B(B — > no charm) ~ 0.02 |14|, |32] is the Standard 
Model branching ratio for charmless decays. Recently, two new measurements of 
the average charm content have been performed. The CLEO Collaboration has 
presented the value n c = 1.16 ± 0.05 |39|, B3[, and the ALEPH Collaboration has 



reported the result n c = 1.23 ± 0.07 [||J]. The average is 

n c = 1.18 ±0.04. (46) 

The naive parton model predicts that Bsl — 15% and n c ~ 1.2; however, it 
has been known for some time that perturbative corrections could change these 
results significantly 0]. With the establishment of the I/tuq expansion, the 
non-perturbative corrections to the parton model could be computed, and their 
effect turned out to be very small. This led Bigi et al. to conclude that values 
£>sl < 12.5% cannot be accommodated by theory, thus giving rise to a puzzle 
referred to as the "baffling semileptonic branching ratio" | 15]| . Recently, Bagan et 



al. have completed the calculation of the 0(a s ) corrections including the effects of 



the charm-quark mass |T6[ , finding that they lower the value of -Bsl significantly. 

The analysis of Bagan et al. has been corrected in an erratum [TO. Here we 
shall present the results of an independent numerical analysis using the same 
theoretical input. As the subject is of considerable importance, we shall explain 
our analysis in detail. The semileptonic branching ratio and n c depend on the 
pole masses of the heavy quarks, which we allow to vary in the range 

m b = (4.8 ± 0.2) GeV , m b - m c = (3.40 ± 0.06) GeV , (47) 

corresponding to 0.25 < m c jm b < 0.33. Here rrib is the pole mass defined to one- 
loop order in perturbation theory. The difference — m c is free of renormalon 
ambiguities and can be determined from spectroscopy (see, e.g., ref. 0). Bagan 
et al. have also considered the theoretical predictions in a scheme where the quark 
masses are renormalized at a scale \i in the MS scheme. We discuss this scheme in 
appendix B. At order 1/ml in the heavy-quark expansion, non-perturbative ef- 
fects are described by the single parameter ^%{B) defined in (|[); the dependence 
on the parameter ^{B) is the same for all inclusive decay rates and cancels out 
in i?sL and n c . Moreover, the results depend on the scale [i used to renormalize 
the coupling constant a s (fi) and the Wilson coefficients c±([i) entering the non- 
leptonic decay rate. They also depend on the value of the QCD scale parameter 
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Aqcdj which we fix taking a s (mz) = 0.117 ± 0.004. The corresponding uncer- 
tainty is smaller than that due to the variation of the mass parameters and is 
added in quadrature. For the two choices fj, — m b and 11 = m b /2, we obtain 



B 



SL 



n r 



12.0 ± 1.0%; 
10.9 ± 1.0%; 

1.20 =F 0.06; 

1.21 =f 0.06; 



V = 

V = 



= m b , 
= m b /2, 
m b , 
m b /2. 



(4* 



The errors in the two quantities are anticorrelated. Notice that the semileptonic 
branching ratio has a much stronger scale dependence than n c . This is illustrated 
in figure |2], where we show the two quantities as a function of fx. By choosing a 
low renormalization scale, values I?sl < 12.5% can easily be accommodated. The 
experimental data prefer a scale \xjm b ~ 0.3-0.5, which is indeed not unnatural. 
Using the BLM scale-setting method E5 



Luke et al. have estimated that fx 



0.3m b is an appropriate scale in this case 



LI 




m c /m b = 0.33 ■ — 
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Figure 2: Scale dependence of the theoretical predictions for the semileptonic 
branching ratio and n c . The bands show the average experimental values. 

The combined theoretical predictions for the semileptonic branching ratio and 
charm counting are shown in figure [| They are compared with the experimental 
results obtained at the T(4s) and at the Z° resonance. It was argued that the 
combination of a low semileptonic branching ratio and a low value of n c would 
constitute a potential problem for the Standard Model [42]. However, with the 
new experimental and theoretical numbers, only for the low-energy measurements 
a small discrepancy remains between theory and experiment. Note that, using 
([45]), our results for n c can be used to obtain predictions for the branching ratio 
B(B — > X C cg/), which is accessible to a direct experimental determination. Our 
prediction of (22 ±6)% for this branching ratio agrees well with the preliminary 
result of the CLEO Collaboration: B(B -> X ccV ) = (23.9 ± 3.8)% 0. 

Having discussed the status of the theoretical predictions obtained to order 
l/m b in the heavy-quark expansion, we now investigate the spectator contribu- 
tions to the semileptonic branching ratio and n c . This extends, in the context 
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*8 9 10 11 12 13 14 
B SL (%) 

Figure 3: Combined theoretical predictions for the semileptonic branching 
ratio and charm counting as a function of the quark-mass ratio m c /m), and 
the renormalization scale /x. The data points show the average experimen- 
tal values for B^l and n c obtained in low-energy (LE) and high-energy (HE) 
measurements. 



of the heavy-quark expansion, the phenomenological study presented in ref. [[LJJ]. 
We consider the average of I?sl and n c for B~ and mesons,^ and write for the 
spectator contributions to these quantities 

A5 S L,s P ec = hBi + b 2 B 2 + 6 3 £i + b 4 e 2 , 
An CjSpec = n x B x + n 2 B 2 + n 3 e 1 + n i e 2 . (49) 

The coefficients 6, and rij are given in table [5|. If, as previously, we eliminate E\ 
from these equations using the constraint ([$3|) imposed by the measurement of 
t(B~) /r(i?d), and we allow that the parameters Bi take values in the range to 
2, we obtain 

ABsL.spec ^ (-2.l£ 2 + 0.2 ±0.3)%, 

An CiSpcc ~ (1.2 ± 0.1) ASgL.spec • (50) 

For reasonable value of e 2 , we expect a contribution to the semileptonic branching 
ratio of order 1% or less, and a negligible effect on n c . However, without a detailed 
calculation of the hadronic parameters we cannot obtain a quantitative prediction 
of the spectator contributions. Nevertheless, we find it interesting that there is 
at least a potential to change, and in particular to lower, the value of B$l by 
0.5-1%. To achieve such a decrease requires that the hadronic parameter e 2 , 
which parametrizes the matrix element of the colour octet-octet operator T|_ P 



3 These are approximately the quantities measured experimentally. However, measurements 
at LEP receive a contamination from B s and 6-baryon decays. 
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in fll3|), is positive and of order 0.3-0.5.Q It will be interesting to see if future 
calculations of this parameter will confirm or rule out this scenario. 

Table 5: Coefficients bi and rii (in %) appearing in (f4"9"|). 





h b 2 b 3 6 4 


m b /2 
m b 

2171b 


0.35 -0.02 -2.60 -1.37 
0.19 -0.03 -2.56 -1.42 
0.03 -0.05 -2.49 -1.42 


fi 


n x n 2 n 3 n 4 


m b /2 
m b 
2m b 


0.50 -0.03 -4.17 -1.54 
0.25 -0.03 -3.79 -1.44 
0.03 -0.05 -3.51 -1.36 



6 Conclusions 

In this paper, we have studied spectator effects in inclusive decays of beauty 
hadrons. Although these effects are suppressed by three powers of Aq CD /m b in 
the heavy-quark expansion, they cannot be neglected because of the large phase- 
space factor for two-body scattering. The contributions of spectator effects to 
inclusive decay rates are given by the hadronic matrix elements of the four local 
operators in (0). For mesons, we have expressed these matrix elements in terms 
of the hadronic parameters B\i and defined in (|i~4"|). For the baryon, heavy- 
quark symmetry reduces the number of independent matrix elements from four 
to two, which we parametrize by r and B as defined in ([27]) and fl28l) . Although 
our parametrization is motivated by commonly made simplifications, such as the 
vacuum insertion and the valence-quark approximations, we stress that it is intro- 
duced without any loss of generality. For a complete understanding of spectator 
effects, it will be necessary to evaluate these parameters non-perturbatively, e.g. 
in lattice simulations. 

We find that in predictions for lifetimes and the semileptonic branching ratio 
of B mesons, the coefficients of the colour octet-octet non-factorizable operators 
are much larger than those for the colour singlet-singlet factorizable operators. 
Thus the contributions from the non-factorizable operators cannot be neglected, 
even though their matrix elements are suppressed in the large- N c limit. 

The ratio T(B~)/r(Bd) is particularly sensitive to non-factorizable contribu- 
tions [see ( |3~2"D and table [3|, making it difficult to predict this quantity with a 

4 We recall, however, that according to ( |39| ) a positive value of £2 increases the theoretical 
prediction for the lifetime ratio T(A b )/T(Bd). 
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precision of better than about 10%. For example, assuming that the magnitudes 
of the parameters e% and £2 are smaller than 0.1 or 0.2, we find that the predic- 
tions for this lifetime ratio lie in the ranges 0.93-1.11 and 0.84-1.20, respectively.^ 
However, in our opinion, even if the experimental result had been outside these 
ranges, the most likely explanation would have been that the parameters are 
larger, rather than a failure of the heavy-quark expansion. The experimental 
measurement of t(B~)/t(B^) imposes the constraint (^) upon the parameters, 
which allows us to eliminate ex in other relations. On the other hand, within the 
heavy-quark expansion there is only room for a very small deviation of the ratio 
t(B s ) /t(Bcl) from unity due to SU(3)-breaking effects. We estimate these effects 
to be of order 1%. 

Understanding the low experimental value of the lifetime ratio r(Ab)/r(fid) 
remains a potential problem for the heavy-quark theory. If the current exper- 
imental value persists, there are two possibilities: either some hadronic matrix 
elements of four-quark operators are significantly larger than naive expectations 
based on large- A^, counting rules and the quark model, or (local) quark-hadron 
duality fails in non-leptonic inclusive decays. In the second case, the explana- 
tion of the puzzle lies beyond the heavy-quark expansion. In the first case, it is 
most likely that the baryonic parameter r is much larger than most expectations 
based on quark-model estimates. It will be interesting to see whether future, 
field-theoretic calculations will yield values of r which are sufficiently large. Until 
such calculations become available, a conventional explanation of the A^-lifetime 
puzzle cannot be excluded. 

Finally, we have performed an analysis of the semileptonic branching ratio 
of the B meson (-Bsl) and of the average number of charmed particles produced 
per decay (n c ). Our results are summarized in figures and [3[ There is a 
significant dependence on the predictions for the semileptonic branching ratio 
on the renormalization scale fi, which is a manifestation of our ignorance of 
higher-order perturbative corrections. The results for n c , on the other hand, 
are almost independent of fi. This scale dependence weakens considerably the 
anticorrelation in the theoretically allowed values for B^h and n c observed in 



ref. P5 |. In our view, given the theoretical uncertainties and the disagreement 
between the experimental values for the semileptonic branching ratio obtained in 
low- and high-energy measurements, there is at present no discrepancy between 
theory and experiment for 1?sl and n c . We have also studied the contributions 
of spectator effects for these quantities and find that they are negligible for n c , 
whereas they can potentially change the prediction for Bsl by up to about 1%. 

Note added: After completing this work we became aware of a paper by 
I.I. Bigi (preprint UND-HEP-96-BIG01, June 1996 ||hep-ph/9606405| ), who dis- 



5 Alternatively, if we assume that the magnitudes of £i(Mhad) an< l e 2(Mhad) renormalized at 
a hadronic scale are less than 0.1 and 0.2, then, using (f42|), we find that the corresponding 
predictions lie in the ranges 0.96-1.14 and 0.87-1.23. 
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cusses theoretical predictions for beauty lifetimes, making strong claims concern- 
ing the theoretical predictions for the lifetime ratio t(B~)/t(B^). In view of our 
discussion in section [OL we must disagree with some statements made in this 



paper. 
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Appendix A: Renormalization-group evolution 



The four-quark operators appearing in the heavy-quark expansion are convention- 
ally renormalized at the scale \x = m b . However, one may use the renormalization- 
group to rewrite them in terms of operators renormalized at a scale /x 7^ mj,. The 
renormalization-group evolution is determined by the anomalous dimensions of 
the four-quark operators in the HQET, where the b quark is treated as static 
quark ||25|| . In the literature, this evolution is sometimes referred to as "hybrid 
renormalization" |TD, E7L |23. 



We find that the operators O v _ A and Ty_ A , and similarly O q s _ P and T|_ 



mix under renormalization. At one- loop order, the mixing within each pair (O, T) 
is governed by the anomalous dimension matrix 



7 



2vr 



C F 
C F 



1 \ 
1 



(A.l) 



V 2N r 2N r J 



which has eigenvalues and 3N C . Here N c is the number of colours, and Cp = 
— 1)/2N C is the eigenvalue of the quadratic Casimir operator in the funda- 
mental representation. The operators defined at the scale m b can be rewritten in 
terms of those defined at a scale \i 7^ m b . In leading logarithmic approximation, 
the result is 



0(m b ) 
T(m b ) 



0{pL)- — {K-\)T{ji) 

T(ji) - §§ (« - 1) 0{jx) , 



(A.2) 



where 



"=(^1 • (A.3) 

\a s {m b ) J 

and p = y N c — |n/ is the first coefficient of the /3-function (n/ = 3 is the 
number of light quark flavours). 

Given the evolution equations ( A.2|) for the four-quark operators, it is imme- 
diate to derive the corresponding equations for the hadronic parameters defined 
in (|TJ), (fZ7|) and We obtain: 



3iV c /2/3o 



Bi(m b ) 
£i(m b ) 



2C p . 

1 + 



1 + 



Bi(^) ~ («- 1) 
C 



r(m 6 ) = « r(/i) + — (k - 1) £(//) r(/i) , 
B(m 6 ) r(m 6 ) = r(/i) . 



(A.4) 
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Of course, introducing parameters renormalized at a scale \i 7^ m b would simply 
amount to a reparametrization of the results and as such is not very illuminating. 
However, the evolution equations are used in section [| to study the sensitivity of 
our results to unknown higher-order corrections. 

As an illustration, we study the evolution from \x = m b down to a typical 
hadronic scale /ihad *C mj, which we choose such that a s (/ihad) = 0.5 (corre- 
sponding to /^had ~ 0.67 GeV). We find 

Bi(m b ) ~ 1.485; (/i had ) - 0.36^ (/i had ) , 

£i(m b ) ~ 1.06£i(Aihad) - 0.085; (/ihad) , 

r(m b ) ~ [1.54 + 0.185(/i h ad)]r(/ihad), 

S(m 6 ) ~ , (A.5) 

1.54 + 0.185(/i had ) 

indicating that renormalization effects can be quite significant. If one assumes 
that the matrix elements renormalized at the scale /ihad can be estimated using 
the vacuum insertion hypothesis for mesons and the valence-quark approximation 
for baryons, then 

fl<(/w) - f 4rr4 = K ~ 8/9 - °- 68 > ^(axw) - , 5(/i had ) ^ 1 , 

(A.6) 

where the factor k 8//9 in the first equation arises from the anomalous dimension 
of the axial current in the HQET [27, Using ( |A.4|) , we then find that the 



parameters defined at a renormalization scale m b (as used throughout this paper) 
would be Bj(m b ) = B 2 {m b ) ~ 1.01, Si{m b ) = e 2 {m b ) ~ -0.05, r(m b )/r(fi had ) ~ 
1.72, and B(m b ) ~ 0.58. Thus, for mesons the violations of the factorization 
approximation induced by the evolution from /ih a d up to m b remain small, i.e. we 
find Bi (777.5) ~ 1 and Si{m b ) ~ 0. We stress, however, that we do not want to 
suggest that the choice of parameters in ( |A.6| ) is actually physical. 



For mesons, the large- N c counting rules jig, f|9| imply that the parameters £j 
are of order l/iV c , whereas the parameters Bi are of order unity. These results 
are respected by the evolution equations (|A.4 ), which in the large- N c limit take 
the form 

B l {m b ) = ^ Bi{fi) + 0{1/ N c ) , 

e,(m b ) = £i (fi) - ^- (^00 - 1) ^(/i) + 0(l/iV c 2 ) , (A.7) 

where /too = [a s (fj>) / a s (m b )] 9 / 22 . Under a change of the renormalization scale, the 
parameters e$ stay of order 1/N C , whereas the parameters Bi change by a factor 
of order unity. 
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Appendix B: £>sl an d n c in the MS scheme 

The semileptonic branching ratio and n c can also be calculated using running 
quark masses renormalized in the Ms scheme rather than pole masses. To compare 
the results in such a scheme to those presented in our work, we have to relate 
the ratio of the pole masses to the ratio of the running masses. There is some 
freedom in how to do this translation. Since in the expressions for the (partial) 
inclusive decay rates radiative corrections are included to order a s (fi) only, it is 
consistent to work with the one-loop relation 

mc = m c {fi) ( _ 2a s {n) ^ m e (/i) \ 
m b m b (n) \ 7T m b (n) ) ' 

We shall refer to this choice as scheme Msl. Alternatively, one may prefer to 
resum the leading and next-to-leading logarithms to this relation, which leads to 

m c (/i) _m c ( a s (m c )y o/2/3 ° ( 1 a s (m c ) - a s (m b ) (^fa - *y Q p x 4\ , ^ 



m b (fi) m b \a s (m b )) I vr V 8/3^ 3 



where 70 and 71 are the one- and two-loop coefficients of the anomalous dimension 
of the running quark mass, and (3q and f3\ are the coefficients of the /^-function. 



We shall call this scheme ms2; it has been adopted in the work of Bagan et al. [[T6 |. 
Our results for these two versions of the Ms scheme are: 

5 SL (MS1) 

n c (M51) 

(B.3) 



w c (ms2) = 

and the combined predictions for Bsh and n c are shown in figure |4j. Contrary 
to the case of the on-shell scheme, the calculations in the Ms scheme become 
unstable for low values of the renormalization scale. For this reason, we only 
present result for fi > m b /2. The results obtained in the scheme Msl are close to 
those obtained in the on-shell scheme and presented in section [|. In the scheme 
Ms2, on the other hand, we find lower values for and higher values for n c . We 
note that our results for this scheme do not coincide with the numbers presented 
in the erratum of ref. [ 16| ; in particular, we do not find the large values of n c 



11.6 ±0.9%; 


/' 


= m b , 


10.7 ±0.9%; 


/' 


= m b /2, 


1.20 =F 0.06; 


n -- 


= m b , 


1.20 =F 0.06; 


n -- 


= mb/2, 


10.9 ±0.9%; 


/' 


= m b , 


10.3 ±0.9%; 


/' 


= m b /2, 


1.25 =F 0.05; 


li -- 


= m b , 


1.24 =f 0.06; 


n -- 


= m b /2, 



reported there. The numerical differences are mainly due to the fact that Bagan 
et al. use lower values for the charm-quark mass (they use m c = 1.33 GeV for 
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the central value of the pole mass rather than 1.4 GeV), and that they multiply 
each partial decay rate by the numerical factor 



1 + 



^-51n^ 
3 Mn ^ 



+ 



(B.4) 



which we omit since it cancels trivially in the dimensionless quantities -Bsl and 
n c . Note that, in particular, this factor is responsible for the large apparent scale 
dependence of the results presented in ref. ||16|| . Another difference is that we 
include in the calculation an estimate of the contribution from charmless decays, 
(2±1) 



using B no charm 
a factor 0.99 and n r 



T|, 0, H and IslVisl c ~ 1%. This lowers 5 S l by 



by a factor 0.98. 



1.4 
1.3 
c°1.2 



1.1 
1 

8 9 10 11 12 13 14 
B SL (%) 

Figure 4: Combined theoretical predictions for the semileptonic branching 
ratio and charm counting as a function of the quark-mass ratio m c /mb and the 
renormalization scale fi. The solid lines refer to the scheme MSl, the dashed 
ones to Ms2. 




It may be argued that the apparent large scheme dependence of the results 
for the semileptonic branching ratio and n c prevent a reliable theoretical predic- 
tion. However, as we have shown above the main reason is that the numerical 
value of the quark mass ratio m c /m(, can be quite different in different schemes 
(m c (/i)/mfe(yu) ~ 0.8m c /mb in the scheme MSl, and 0.7 m c /mb in Ms2). Since 
the dependence of -Bsl and n c on the quark-mass ratio comes simply from phase 
space (and is particularly strong for the channel h — > ccs), we feel that the on- 
shell scheme is more adequate for performing the calculation. In other words, we 
expect that in the MS scheme one would encounter larger higher-order corrections, 
once the calculation is pushed to order a 2 s and higher. 
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